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Abstract 



We derive the 1-point probability density function of the smoothed Abell-ACO cluster density field 
and we compare it with that of artificial cluster samples, generated as high peaks of a Gaussian 
field in such a way that they reproduce the low-order (2- and 3-point) correlation functions and 
the observed cluster selection functions. We find that both real and simulated pdfs are well ap- 
proximated by a log-normal distribution even when the Gaussian smoothing radius is as large as 
40 h -1 Mpc. Furthermore the low-order moments of the pdf are found to obey a relation 7 oc a 4 , 
with 7 being the skewness. Since clusters have not had enough time to significantly depart from 
their original birth-place positions, these results are consistent with them being high-peaks of an 
underlying initial Gaussian density field. 

A by-product of our analysis is that when we rescale the pdf cluster moments to those of the 
QDOT-IRAS galaxies, using linear biasing with b c j ~ 4.5 and for the common smoothing radius of 
20 h~ l Mpc, we find them to be significantly smaller than those directly estimated from the QDOT 
data by Saunders et al. (1991). 
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1 Introduction 



The study of the large-scale galaxy distribution is of fundamental importance in modern cosmol- 
ogy. A complete knowledge of the present cosmic density fluctuation field 5 m = , in the linear 
regime, would yield informations on the shape of the spectrum of primordial perturbations which 
gave rise, in the gravitational instability picture, to the present cosmic structures and their distri- 
bution. Within the framework of the most common galaxy formation theories, the standard two 
assumptions are (i) that the initial 5 m field is described by a Gaussian distribution and (ii) that 
some sort of biasing mechanism is at work. In the most popular biasing models, galaxies form at 
the high peaks of the underlying matter field (Kaiser 1984; Bardeen et al. 1986). This model was 
introduced to reconcile, among other things, the observed low values of £1 with the inflationary 
model (Olive 1990). 

Many statistical approaches have been used to study the distribution of extragalactic structures. 
The most common one is the two-point correlation function, £(r), (cf. Peebles 1980) which has 
been found to have the form: 

f (r) = (r/r )~\ (1) 

with r D ss 5 h~ l Mpc and 7 ss 1.8 from 10 h' 1 Kpc up to 10 h' 1 Mpc (Davis & Peebles 1983). 
Cluster of galaxies exhibit similar kind of correlations up to ~ 50 h~ l Mpc with the same slope and 
r D ss 20 - 25 h' 1 Mpc (cf. Bahcall & Soneira 1983; Klypin & Kopylov 1983). Another statistical 
measure is the 3-point correlation function, for which a non-zero signal implies non-Gaussian 
statistics. The analysis of both 2-D and 3-D samples has shown that galaxies do have a non-zero 
connected 3-point correlation function with a hierarchical coefficient, Q ~ 1 (Groth & Peebles 1977; 
Sharp, Bonometto & Lucchin 1984). Similar results hold also for clusters of galaxies (Toth et al. 
1989, Jing & Zhang 1989, Jing & Valdarnini 1991 [hereafter JV]; Borgani, Jing k Plionis 1992). 

These results are of great interest, since non-Gaussianity might arise either from gravitational 
clustering, or from biasing mechanisms, or finally from intrinsic non-Gaussian statistics of the ini- 
tial matter field. A complete knowledge of the N-point correlation function would allow one, using 
the linear biasing assumption, to recover the statistical properties of the S m field. For N > 3, how- 
ever, this approach becomes impractical due to the limited extension of the available extragalactic 
samples. 

An alternative method, which is becoming increasingly popular, is the study of the probability 
density function (pdf). Theoretically the knowledge of the N-point correlation function is equivalent 
to that of the pdf and its moments, but for real data sets the latter are easier to measure. Because 
of the discrete nature of the available samples of extragalactic objects, and in order to obtain a 
continuous density field one is forced to smooth the discrete distribution with some function (usually 
a top-hat or a Gaussian one is used). For a normally distributed primordial density field and as 
long as the variance a 1 = (5 2 ) is small, under the action of gravity, the deviation of the initial pdf 
from the Gaussian shape is well approximated by the Edgeworth expansion (Colombi 1994). 

In the non-linear regime N-body simulations have shown the the pdf becomes non-Gaussian as 
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the clustering evolves (Bouchet, Schaffer & Davis 1991; Bouchet & Hernquist 1992; Lahav et al. 
1993; Juszkiewicz et al. 1993; Kofman et al. 1994). For a <C 1 the shape of the pdf is well fitted 
by a lognormal distribution. The pdf of different galaxy samples has been estimated by Bouchet 
et al. (1993), Gatzahaga and Yokoyama (1993). Kofman et al. (1994) have compared the pdf 
derived form CDM N-body simulations with that of the IRAS sample and the one recovered using 
the POTENT procedure with = 1. Their main conclusion is that, if galaxies trace the mass, the 
observed pdf can be fitted using Gaussian initial conditions. 

Furthermore, the time evolution of the moments of the pdf has been studied both analytically 
and with N-body simulations (Goroff et al. 1986; Coles & Frenk 1991; Juszkiewicz et al. 1993; 
Juszkiewicz, Bouchet & Colombi 1993; Kofman et al. 1994; Colombi 1994). For a given filter and 
smoothing radius, the normalized skewness of the matter distribution, S3 = (S 3 }/(5 2 ) 2 , will depend 
on the shape of the power spectrum P(k). Juszkiewicz , Bouchet & Colombi (1993), using second 
order theory, found S3 ~ 34/7 — (3 + n) for a top-hat filter and scale-free spectra. These results 
agree with those obtained from N-body simulations (Kofman et al. 1994). The proportionality 
(5 3 ) oc o holds also at a fixed time when one changes the smoothing radius (Coles &; Frenk 1991; 
Kofman et al. 1994). 

In this paper we estimate the pdf for the distribution of Abell+ACO clusters of galaxies (Abell 
1958; Abell, Corwin & Olowin 1989). Since clusters of galaxies are the brightest objects in the 
sky they can probe the very large scales, reliably up to ps 300 h^ 1 Mpc. On these scales the 
density fluctuations are well within the linear regime and the deviations of the pdf from a Gaussian 
distribution is determined either by the fact that they are biased tracers of the underlying field 
or due to non-Gaussian initial conditions. The main aim of this paper is to establish whether the 
present cluster pdf can be derived entirely under the assumption that clusters form at the high 
peaks of a Gaussian background field. The basic procedure is to compare the observed cluster pdf 
with that measured from the distribution of simulated clusters, generated as peaks of a initially 
Gaussian background field and having the same boundaries and selection effects as the real data. 

2 Data Samples 

2.1 Real cluster samples 

We apply our analysis to the combined Abell-ACO R > cluster sample, as defined in Plionis 
& Valdarnini (1991) [hereafter PV91] and analysed in Plionis, Valdarnini & Jing (1992) [hereafter 
PVJ]. The northern sample, with dec> —17° (Abell), is defined by those clusters that have measured 
z < 0.1, while the southern sample, (ACO) with dec< —17°, is defined by those clusters with 
"Ho < 16.4 (note that with this definition and due to the availability of many new cluster redshifts 
only 7 ACO clusters have mio estimated redshifts from the mio — z relation derived in PV91). Both 
samples are limited in Galactic latitude by \b\ > 30°. The redshifts we have used are mainly from 
Struble & Rood (1987) and Postman, Huchra & Geller (1992) with additions from taken Rhee & 
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Katgert (1988) and Batuski et al. (1991), the original ACO paper, Vettolani et al. (1989); Cappi 
et al. (1991) and Muriel et al. (1991) as well as by our own cross-correlation analysis of the Fairall 
& Jones (1991) galaxy redshift catalogue with the cluster catalogues. The total number of clusters 
in our samples is 357 Abell and 157 ACO ones. 

We take into account the effect of Galactic absorption, which is assumed to follow the usual 
cosecant law: 

P(|6|) =dex [a(l-csc|6|)] (2) 

with a ~ 0.3 for the Abell sample (Bahcall & Soneira 1983; Postman et al. 1989) and a ~ 0.2 for 
the ACO sample (Batuski et al. 1989), by weighting each cluster by 1/P(|6|), a correction which 
explicitly assumes that the unobserved (obscured) clusters are correlated with the observed ones, 
a reasonable assumption under the well known correlation properties of clusters. 

We have derived the cluster-redshift selection function, P(z), by fitting the cluster density, as 
a function of z (cf. Postman et al. 1989; PVJ), by: 

m\ =1 . . . f for \- 2c (3) 

[ oc exp(— z/z ) for z > z c 

where z c is the maximum redshift at which the sample exhibits a roughly constant density (for more 
details see PVJ). We find z c « 0.079 and 0.066 for the Abell and ACO samples, which correspond 
to R ~ 230 and 190 h~ l Mpc respectively. Cluster comoving distances are estimated using the 
standard relation (Mattig 1958): 

R = H q{{l+z) ( q ° Z q ° ){1 ~ V / 2^+l)) W 

with H = 100 h km sec" 1 Mpc" 1 and q Q = 0.2. 

The Abell and ACO cluster densities, out to their limit of completeness (z c ), is ~ 1.4 x 10" 5 h 3 
Mpc -3 and ~ 2.1 x 10 -5 h? Mpc -3 , corresponding to a characteristic length scale of ~ 41 

h~ x Mpc and 36 h~ x . The higher space-density of ACO clusters is partly due to the huge Shapley 
concentration (Shapley 1930; Scaramella et al. 1989; Raychaudhury 1989), but a significant part 
is also due to systematic density differences between the Abell and ACO cluster samples, as a 
function also of z, which have been noted in a number of studies (cf. PV91 and references therein) 
and which could be attributed to the high sensitivity of the Illa-J emulsion plates. 

Since we want to use a cluster sample covering the whole sky we need to take into account the 
density variations of the Abell and ACO cluster samples. We will do this by renormalizing the ACO 
density to the Abell one using a radial matching function, W(R), as in PV91. Since in this paper 
we are mostly interested in the cluster density fluctuations Sp/p, it is equivalent to normalize the 
density of the Abell to the ACO or vice-versa. In order to test the robustness of the results in the 
choice of the W(R) function, we will use a number of such functions, among which the W(R) = 1 
(no matching) and the W(R) = constant (no radial dependance) cases. 

Note that in PVJ we had estimated the 2-point correlation function for our samples of clusters 
and we found that the slope of the derived 2-p function has a value ~ 1.8 ± 0.2 for both Abell 
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and ACO samples. However their amplitudes are slightly different, with Abell clusters having 
r Q ~ 18 ± 4 Yr x Mpc (bootstrap errors used) out to < 50 h" 1 Mpc while the ACO clusters have 
r Q ~ 22 ± 10 Mpc but only out to ~ 30 h~ x Mpc. This could be, however, due to the fact that 
the ACO sample is defined in a relatively small solid angle and therefore the large wavelengths may 
be undersampled. Based on our previous analysis (Jing et al. 1992), we believe that the cluster 
correlations are reliable and mostly unaffected by cluster contamination effects (c/. Sutherland 
1988). 

In PVJ we also estimated the 3-point correlation function using the moment method (cf. Sharp, 
Bonometto Sz Lucchin 1984; Jing & Valdarnini 1991; Borgani, Jing & Plionis 1991). We found a 
good fit to the spatial 3-point function provided by the hierarchical expression: 

C(n,r 2 ,r 12 ) = Q[e(r 1 )^(r 2 )+e(ri)e(r 12 )+e(r 2 )^(r 12 )], (5) 

with Q Abell ~ 0.6 and Qaco ~ 0.7, with relative bootstrap error of ~ ±0.2. 

Finally, in order to minimize the uncertainties due to the approximate character of the redshift 
selection function, P(z), and of the radial function, W(R), especially at large distances, we will 
restrict our analysis to R max = 240 h~ x Mpc. 

2.2 Simulated cluster samples 

The procedure used to generate our artificial cluster catalogues has been presented in PVJ. Here 
we only remind the reader the main steps of our procedure. 

We generate our simulated cluster catalogues using a method similar to that of Postman et al. 
(1989). We identify the position x of points corresponding to clusters of galaxies according to the 
following prescription: N p points are randomly placed in a cube of size 640 h- 1 Mpc with N g = 64 3 
grid points. We keep points only within the sphere of radius r max = 320 h^ 1 Mpc centered at the 
cube centre. The points are then displaced from their original positions x r using appropriately the 
Zel'dovich approximation. The power spectrum of the density fluctuations, P(k), has a Gaussian 
distribution with random phases. P(k) is chosen so that the perturbed particle positions give a 
correlation function in the desired range. In accordance with Postman et al. (1989), we use the 
following power spectrum: 

P(k) = Ak n exp(-|fc| 2 /A 2 ) 9(|jfe|) (6) 

where A -1 = 0.1 Yr 1 Mpc, 6(|fe|) = for |fe| > 27r/80 h~ x Mpc, A is a normalization constant and 
n is the spectral index. 

We assign 'clusters' to the peaks of the background field according to the following prescription: 
each particle inside the cube is assigned a value v such that 5 S > ua, where a is the rms density 
fluctuation within the cube and g is the nearest grid point to the particle. Here 5 S is the density 
field smoothed with a Gaussian filter function having a smoothing radius of 10 h^ 1 Mpc. The 
parameters A and n and v must be chosen such that the peak-peak correlation function in the 
range 10 - 60 h~ l Mpc, corresponds to the observed £ cc . We find that for A m 1.9 x 10 6 , n« —1.5 
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and v ~ 1.3 we obtain a distribution of simulated clusters with {r ) s im = 21 h^ 1 Mpc. Furthermore 
these clusters have by construction non-zero higher order correlations and we have found in PVJ, 
using the direct triplets-counting method, Q s i m — 0.6 ±0.16, where the uncertainty is the ensemble 
one. 

Finally we note that our simulated cluster catalogues have the same geometrical boundaries, 
the same redshift and Galactic obscuration selection functions, the same low-order correlation 
properties as well as the same mean space density as the real clusters. 



3 Method 

In order to obtain a continuous cluster density field we smooth the cluster distribution in a 24 3 
cube (480 3 /i" 3 Mpc 3 ) using a Gaussian kernel: 

W(x t - x g ) = (27TUL) ~ 3/2 exp (-^^) (7) 
The smoothed cluster density, at the grid-cell positions then: 

P{X * } - /W(|*-* g |)d 3 z (8) 

where the sum is over the distribution of clusters at positions Xj. In order to study the cluster 
probability density function, P(p), at different scales, we use four smoothing radii; R sm = 20, 30, 
40 & 50 h- 1 Mpc with |x; - : Kg I < 3i? sm . In this case the integral at the denominator of eq.(8) is 
~ 0.97. There are, however, two problems which we have to resolve: 

• Due to the geometrical boundaries of the observation volume and the zone of avoidance, each 
Gaussian sphere is incomplete. Using a Monte-Carlo method, we estimate this incompleteness, 
by defining a grid completeness factor. 



J n(x)p(x)W(\x — x g |)d 3 x 
/ n(x)W(|x - x g |)d 3 x 



where n(x) is the density of a random distribution of points (thus n(x) = constant) and 
p(x) represents the selection functions and geometric boundaries. Note that by definition 
< f(x s ) < 1. The smoothed cluster density is then corrected, at each grid point x s , by 
l//(x g ). However, to reduce the uncertainty due to the approximate nature of this correction 
we use in our analysis only those grid-cells with f(x s ) > 0.8. Due to the small number of 
such cells for large smoothing radii, we use also f(x g ) > 0.7 but only for the R sm > 40 cases. 

• Discreteness effects can be introduced for small R sm , if the number of cluster counts, in the 
Gaussian sphere, is small or if the smoothing fails to create a continuous density field. The 
latter occurs only for R sm = 20 h^ 1 Mpc where some cells have p sm = 0. This biases the P{p) 
at p/(p) <C 1. This brings up the issue of which is the optimal smoothing radius to derive 
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the P(p) of a discrete distribution of points. From considerations presented in the following 
section we find that the optimal radius is of the order the mean interparticle separation, ie. 
~ 30 - 40 h- 1 Mpc for the Abell+ACO clusters. 

By applying the selection functions and the geometric boundaries of the real cluster distribution 
to that of the simulated clusters we have tested the robustness of our correction procedure by 
comparing bin- by-bin the original simulation cluster density contrast, S or i (whole cube), with the 
recovered density contrast, 5 rec . We find an excellent agreement (Figure 1) and we have also 
verified that the pdfs, before and after the application of the selection functions and the correction 
procedure, are identical. 



4 Results & Discussion 



In figure 2 we present the simulation-cluster mean pdfs (over the 50 simulations), for all four 
smoothing radii, after the application of the selection functions and the subsequent corrections, 
discussed in the previous section. The errorbars are the la scatter around the ensemble mean. We 
plot the Gaussian fit (dashed line) to the cluster pdf, given by: 
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where a is the standard deviation of g(= p/(p)), and we also plot the log-normal fit, which has 
been found to give an extremely good fit to the CDM density and the IRAS galaxy pdf in the 
weakly- linear regime (Kofman et al. 1994). The log- normal distribution is given by: 
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(11) 



where pi and is the mean and standard deviation of ln£. 

From figure 2 it is evident that the log- normal distribution fits well the simulation cluster pdf for 
all four R sm , while the Gaussian one is rejected and although it becomes increasingly a better fit to 
the cluster P(g), as R sm increases, it always provides a worse fit than the corresponding log-normal 
distribution, even for R sm = 50 h^ 1 Mpc. Note, however, that for R sm = 20 h^ 1 Mpc and g < 1 
the log-normal distribution does not fit well the pdf. The reason is that at this smoothing scale and 
range of g there are significant shot-noise contributions (see discussion in previous section). 

The log-normal distribution has been argued (Coles & Jones 1991) to describe the distribution 
of density perturbations resulting from Gaussian initial conditions in the weakly non-linear regime. 
It has been shown to fit well the Lick Euler-Poincare characteristic (a measure related to the genus) 
and the 2-dimensional Lick counts in cells at a resolution which corresponds, at the characteristic 
depth of the Lick catalogue, to scales of ~ 4 — 10 h^ 1 Mpc (Coles & Plionis 1991). However, 
Bernardeau & Kofman (1994) have shown that the log-normal distribution is not a universal form 
of the cosmic density pdf 'but a very convenient fit only in some portion of the (<7,n)-plane (ie, a <C 1 
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and n~ —1)- In this context it seems that it is not a coincidence that the simulation cluster pdf is 
fitted quite well by a log-normal distribution but it could be a consequence of the fact that we used 
a power spectrum with exponent near to —1 (n s j m ~ —1.5). Alternatively, it could be a generic 
feature of high-peak biasing. 

In figure 3 we present the one-point pdf of the Abell+ACO smoothed cluster distribution, for 
three R sm values^ and for two different values of f{x s ), the latter to check the effect of the limited 
number of cells used to derive the pdf. The uncertainty here is given by Poisson sampling errors and 
thus they should be considered only a lower limit to the intrinsic uncertainty. A more reasonable 
indication of the scatter is the ensemble errors seen in figure 2. Again we also plot the best Gaussian 
(dashed line) and log-normal (solid line) fits. We verify that also the real cluster smoothed pdf is 
fitted quite well by a log-normal distribution, especially for R sm > 30 h^ 1 Mpc where discreteness 
effects are minimized. 

We now proceed to estimate the moments of the simulated and real cluster distributions, which 
are defined by: 



where 8 = g — 1. Note that if f(8) is a Gaussian then it is defined in an infinite interval and thus 
it always has /(< —1) ^ 0. However by definition 8 > — 1 and therefore f(8) can be a Gaussian 
only in the limit a — ► 0. The gravitational evolution of the 8 field acts in a way to increase the 
variance, a 2 = (8 2 ), and thus in order for 8 > — 1 to hold, f(8) becomes skewed. Note that the 
skewness, 7 = (8 s ) of a Gaussian distribution (defined in an infinite interval) is zero. In the case of 
the 8 field, however, we have that: 



and thus 7 / (unless a — ► 0). We have already introduced the variance and the skewness of 
f{8). In a similar fashion we can define higher order moments; for example the 4 th moment is the 
kurtosis, K = (8 A ) (cf. Lahav et al. 1993; Kofman et al. 1994). 

An important problem in determining the moments of f(8) is the contribution of discreteness 
effects which can dominate in low-density discrete distributions. If the point-like distribution repre- 
sents a Poisson sampling of an underlying continuous density field then the shot-noise contributions 
can be easily corrected (cf. Peebles 1980). However, the galaxy-cluster distribution can be hardly 
considered a poisson sampling since it rather represents a strongly biased sampling of an under- 
lying galaxy density field. Therefore the poisson shot-noise corrections should not be expected to 
provide a reasonable correction of the discreteness effects. In fact, Borgani et al. (1994) found 
for their 2-d cluster analysis that indeed the poisson shot-noise correction resulted in extremely 
noisy moments and they preferred not to correct for shot-noise effects, which resulted in stable, 
well-behaved moments. 

1 We have not used the R srn = 50 Mpc because there are only a few cells with / > 0.7. In the simulation case 
we overcame this problem due to the large number of available realizations. 






(13) 
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Furthermore, the smoothing process itself suppresses the shot-noise effects considerably (cf. 
Gaztahaga & Yokoyama 1993). From eq.(8) it is easy to show that the power spectrum of the 
smoothed cluster density is related to that of the cluster P c i(k) by 

P s (k) = e" fc2 ^ m 

where (I) is the average cluster separation. Then the best representation of P c i(k) is for R sm ~ (I), 
the original field being oversmoothed for R sm 3> (I) and undersampled for R sm <C (I). 

In Table 1 we present the second, third and fourth moment of the cluster pdf as a function 
of smoothing scale, for the simulations as well as the real data. It is apparent that the cluster 
moments are significantly smaller than those of the simulation clusters. Although the latter are 
constructed in such a way to reproduce the correct slope and roughly the correct amplitude of the 
2-point spatial correlation function of the Abell and ACO (R > 0) clusters (see PVJ), there are 
differences in the amplitude. For example the correlation length of the Abell and ACO clusters 
is r a ~ 18 and 22 h^ 1 Mpc, respectively while of the simulated clusters {r ) s i m — 21 h~ l Mpc 
and since the Abell sample is the dominant one we expect that the overall cluster variance will 
be smaller than that of the simulated clusters. In the same direction works also the fact that the 
simulated clusters have a non-zero £(r) at separations larger than those of the real clusters and 
therefore since the smoothing process mixes a lot of different scales, the non-zero £(r) at large r 
will tend to increase the variance. 

As we discussed in the Introduction many authors found, on the basis of a variety of clustering 
scenarios from Gaussian initial conditions, that the skewness is simply related to the variance of 
the mass distribution by: 

l = B(a 2 ) A (15) 

This result holds for a < 1 (cf. Fig. 5 of Kofman et al. 1994) while a changes either due to 
the clustering evolution or due to different smoothing radii. We therefore plot in figure 4(a) the 
skewness, 7, versus the variance, a 1 for the real cluster pdf and in panel (b) for the simulated 
cluster pdf. For the simulated cluster case we plot the values for each of the 50 simulations and for 
each simulation at 4 different smoothing scales (as discussed in the previous section). Similarly for 
the real clusters we plot the values corresponding to different Abell-ACO normalization functions, 
W(R), and for 3 smoothing scales. It is evident that the plotted points in both, simulation and real 
cluster, cases lie on the same line. To check whether eq.( |l5|) is fulfilled we have used a weighted 
least-square fit to determine the slope and the intercept of the line: 

log 7 = A log a 2 + log B (16) 

and we have found for the cluster case: 

A = 2.09 ± 0.05 5 = 1.94 ±0.15, 



N 



4vr r 1 j., x 9 , 



(14) 



9 



and for the real cluster case: 



,4 = 2.04 ±0.06 5 = 1.8 ±0.2, 



which imply that indeed eq.( 15) is fulfilled. This result is also in agreement with the hierarchical 
relation (7 = 3Qa 4 ) with Q ~ 0.6, consistent with previous cluster correlation function analysis 
(cf. JV, PVJ). The parameter B is also called normalized skewness, S3. In general, the normalized 
moments are given by: 

s "-^h < 17) 

where (S n ) c are the cumulants. The normalized skewness, 53, plotted in Figure 5, for both the 
simulation and real cluster pdf has roughly the same value, and is constant as a function of R sm - 
Such a behaviour is consistent with density distributions evolving from Gaussian initial conditions 
(Coles & Frenk 1991; Kofman et al. 1994). Furthermore, eq.( |l5|) is valid also in the high-peak 
biassing model in CDM universes (Coles & Frenk 1991). This is worth stressing since what we 
measure here is not the skewness induced on the density field by gravity but that arising from a 
threshold peak-selection. 

Therefore, the above results, together with the fact that both the simulated and real cluster pdf 
are fitted by a log-normal distribution, imply: 

• the Abell±ACO clusters are compatible with being high-peaks of a Gaussian background and 

• the spectral index of density fluctuation spectrum, from which the cluster distribution has 
emerged, is compatible with ~ —1.5. 

We emphasize that they are 'only ' consistent because a thorough analysis of the expected behaviour 
of peaks of a non-Gaussian background in 3-D has not yet been investigated thoroughly and there- 
fore we cannot exclude the possibility that such peaks may have a similar behaviour. However, 
there are indications from the study of the moments of simulated galaxy distributions emerging 
from non-Gaussian initial conditions (Weinberg & Cole 1992; Coles et al. 1993), that they scale 
in a manner similar to that of distributions emerging from Gaussian initial conditions. In fact, 
Coles et al. (1993) find that in their initially positive skewed models, S3 is also a constant but 
has a higher value than that of the Gaussian models. Initially negative skewed models are difficult 
to distinguish from Gaussian models on the basis of their normalized moments since they produce 
similar values of S3 (see also Table 2 of Weinberg k, Cole 1992). 

A similar analysis for clusters of galaxies has been performed, but only in 2-dimensions (Borgani 
et al. 1994) and it was found that not only the skewness and the variance, in both Gaussian and 
non-Gaussian models, scale in a similar manner but also S3 ~ 1.8 for all the models. Since however 
the projection from 3 to 2 dimensions Gaussianizes the distribution, as expected from the central 
limit theorem, these results may not be representative of the 3-D case. 

It would be interesting to compare the cluster pdf moments with those of galaxies for similar 
smoothing scales. In fact, we have one smoothing scale (R sm = 20 h^ 1 Mpc) in common with the 
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analysis by Saunders et al. (1991) of the QDOT IRAS redshift sample, and it is this smoothing 
scale for which they found large values of the variance and the skewness and claimed that the CDM 
model is ruled out at a 97% confidence limit. Assuming that the linear biassing model is correct, 
on the relevant scales, and using a biassing factor between clusters and IRAS galaxies of b c i ~ 4.5 
(Peacock & Dodds 1994) we predict 

oijred. = K,Wc = 0-0245 ± 0.005 
li, P red. = b-fr c = 0.0049 ± 0.0018 

while the corresponding measured values of the QDOT sample and of the CDM QDOT look-alikes 
(from Saunders et al. 1991) are (a 2 ,j) = (0.0669,0.025) and (0.0192,-0.001) respectively. It is 
evident that the predicted IRAS moments from the cluster moments and from linear biassing are 
significantly smaller than the Saunders et al. (1991) direct estimation but still larger, although 
marginally so, than the CDM ones. Alternatively we can ask whether there is a unique biasing 
value, b c i, for which the rescaled cluster moments, to the IRAS level, are equal to the QDOT-IRAS 
ones. In fact there is; using b c j 2.7 we succeeded in the above task. However, such a low value of 
b c j does not seem to be supported by many different studies (cf. Jing & Valdarnini 1993; Peacock 
& Dodds 1994). 

These results may be viewed as suggestive that the linear biassing model cannot be applied 
straightforwardly to the moments of the pdf. However, Juszkiewicz et al. (1993) have found that 
eq.( |l5|) still holds for a biased population, if the linear fluctuations are Gaussian and the biased 
density is a local function of the matter density (see also Coles & Frenk 1991). 

Jing & Valdarnini (1993) have found that the quadrupole anisotropy, observed in the COBE 
experiment, can be reproduced by the measured Abell cluster and IRAS galaxies power spectra, if 
their bias factors are b c j ~ 3 — 4 and bj ~ 1.5, respectively (see also Peacock and Dodds 1994). 
These results, together with the value S3 ~ 1.8 which we obtain in this work for Abell+ACO clusters 
and using the linear biasing framework, we estimate that for the matter distribution S^m ~ 8 — 12. 
These values look unrealistically large since Goroff et al. (1986) obtain S% m ~ 3 for a CDM spectra 
and a Gaussian smoothing with R sm > 20 Mpc. 

Although, in the range of interest the cluster power spectrum has an effective slope steeper than 
that of the CDM, the values of Ss m cannot exceed an upper limit of ~ 5 obtained from perturbation 
theory (cf. Fig.l of Juszkiewicz, Bouchet & Colombi 1993). Taken together all these results point 
toward a more cautious approach for the linear bias assumption. On the other hand, however, the 
optical cluster dipole (PV91; Scaramella et al. 1991) is found to be well aligned with that of the 
CBR, which means that clusters do effectively trace the underlying gravitational field. 

We therefore suggest that the linear bias assumption is a good approximation for spatial regions 
where 5 > but breaks down in underdense regions where 6 < 0. The skewness is a measure of the 
asymmetry of the distribution with respect the mean and therefore one cannot simply relate the 
moments of the biased field to that of matter field by a constant factor. 
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5 Conclusions 



We have derived the Abell+ACO cluster probability density function at different smoothing scales 
and we found it to be consistent with a log-normal distribution at all scales considered. Similar 
results were obtained from the analysis of synthetic clusters generated as high-peaks of a Gaussian 
background density field having the same low-order correlation properties and selection functions 
as the real clusters. Furthermore we determined the moments of the pdf and we found them 
to scale according to: 7 = S3CX 4 . Such a scaling is expected within the framework of Gaussian 
initial conditions and gravitational instability. However, since the clusters probe the linear field, 
the scaling is attributed to the the fact that they are biased tracers of the underlying density field. 
Furthermore, this scaling is in agreement with the hierarchical clustering model with S3 = 3Q ~ 0.6, 
consistent with previous correlation function analysis (cf. VJ, PVJ). These results, together with 
the arguments of Kofman et al. (1994), indicate that the observed cluster distribution is consistent 
with it being generated from a initial Gaussian fluctuation spectrum. 
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Table 1: The low-order moments of the Abell+ACO (subfix A) and simulated cluster (subfix sim) 
pdfs, as a function of smoothing scale R sm . The uncertainty in the real cluster moments is given 
from the scatter induced by using different W(R) functions (see text), while in the simulation case 
the uncertainties are due to the scatter around the ensemble mean values. Note that we also give 
the values for the fourth moment, K, although in the text we do not discuss it due to the large 
uncertainty associated with its estimation. Note, however, that for the R sm = 20 h~ l Mpc case, 
where we have a relatively more reliable estimation, S^(= ) ~ 5.1 and 5.3 for the real and 

simulated cluster data, respectively. 



Rsm sim Ksim K-sim 

20 0.496±0.057 0.4456±0.154 1.363±0.7 0.85±0.19 1.44±0.68 5.4±3.1 

30 0.136±0.023 0.0327±0.0168 0.365±0.566 0.36±0.11 0.29±0.2 0.67±0.5 

40 0.048±0.013 0.0040±0.005 0.2±0.078 0.095±0.087 0.173±0.154 
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Figure Captions 



Figure 1. 

Comparison, cell-by-cell, of the original simulated cluster fluctuation field, 8 or i, with the one recov- 
ered after applying the real cluster selection functions and geometrical boundaries and using the 
method described in section 3 (for f(x s ) > 0.8). A random subsample of 10%, from each of the 50 
simulations, has been plotted. 
Figure 2. 

The simulated cluster pdf at the 4 smoothing scales, indicated in each panel (filled symbols). The 
errorbars are the scatter around the ensemble mean (over 50 simulations). The solid line is the 
log- normal fit to the data while the dashed line is the Gaussian fit. 
Figure 3. 

The real cluster pdf at the 3 smoothing scales indicated in each panel. Filled symbols correspond 
to /(x g ) > 0.8 while open symbols to f(x s ) > 0.9, for the R sm = 20 and 30 h~ l Mpc cases, 
respectively. For R sm = 40 h^ 1 Mpc they correspond to 0.7 and 0.6, respectively. The errorbars 
are Poisson errors and the line fits are as in figure 2. 
Figure 4. 

Variance versus skewness plot. The open circles, solid squares, open stars and open squares corre- 
spond to R S m = 20, 30, 40 and 50 Mpc respectively. The solid line is the best fitted line given 
by eq.( |l6|). (a) Real clusters (b) Simulated clusters. 
Figure 5. 

The normalized skewness, S3, versus smoothing scale for both real and simulated clusters. 
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